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We study the five-quark system uudds in tiie standard non-relativistic quark model by 
solving tlie scattering problem. Using the Gaussian Expansion Method (GEM), we perform 
the almost precise multi-quark calculations by treating a very large five-body modelspace 
including the NK scattering channel explicitly. Although a lot of pseudostates (discretized 
' continuum states) with ,P = and J'^ = |^ are obtained within the bound-state approx- 

imation, all the states in 1.4 — 1.85 GeV in mass around 0'''(154O) melt into non-resonant 
continuum states through the coupling with the NK scattering state in the realistic case, i.e., 
i-^ ' there is no five-quark resonance below 1.85GeV. Instead, we predict a five-quark resonance 

^! state of J" = i~ with the mass of about 1.9GeV and the width of T ~ 2.68MeV. Similar 

^ l' calculation is done for the four-quark system ccqq {q = u, d) in connection with X''(3872). 

§1. Introduction: importance of the wave function of multi-quarks 

J> ' From the viewpoint of Elementary Particle Physics, the lattice QCD calculation 

■ is one of the most standard approaches to investigate hadrons including multi-quarks, 
i However, it is rather difficult to extract the "wave function" of hadrons in the lattice 

'sj" ' QCD calculation because it is based on the path-integral, where all contributions of 

■ possible states are summed up and only vacuum expectation values can be obtained, 
i Needless to say, the wave function is one of the most important quantities in 

quantum physics, and, of course, also in Quark-Hadron Physics. In particular for 
multi-quark systems, the analysis with the quark wave function is necessary to clarify 
whether the multi-quark hadron is an exotic resonance state or a two-hadron scatter- 
^ ' ing state. Thus, to extract the state information (the wave function) of hadrons, we 

■ need a reliable calculational method for multi-quark systems instead of lattice QCD. 
One of the attractive methods is the constituent quark-model calculation, since it 
seems workable up to several hundred MeV excitation, and its precise calculation for 
multi-quark systems can be done with the Gaussian Expansion Method (GEM).!!''''^ 

In Nuclear Physics, a precise calculational method for bound and scattering 
states of various few-body systems using GEM^'' has been developed by two of the 
present authors (E.H. and M.K.) and their collaborators, and has successfully been 
applied to light nuclei, light hypernuclei, exotic atoms/molecules and so on.13 

Using GEM within the framework of a constituent quark model, Hiyama et al^ 
investigated for the first time scattering and resonance states of the five-quark system 
uudds under the explicit NK scattering boundary condition. It was made clear that 
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there appears no resonance state around the reported energy of 0^(1540).'^'®''^ 
In this paper, we perform the similar precise calculation of the five-quark system 
uudds with J'^ = ^ and J'^ = | using the linear confinement potential, which 
is more appropriate as indicated by lattice QCD calculations^'''^' Furthermore, we 
apply the same method to the four-quark system ccqq in connection with X°(3872j3 
and discuss about possibility of any bound/resonance state near the D'^D*'^ threshold. 

§2. Five-body quark-model calculation for ©+(1540) 

We study the five-quark system uudds by solving the five-body Schrodinger 
equation {H — E) ^jttm = including the NK scattering channel explicitly, as was 
done in Ref.2). We here take a standard non-relativistic quark- model Hamiltonian 
with the linear- type confining potential V^onf ^-nd the color-magnetic potential Vcm^ 

V,oAn,) = ^(^a n, + vo), VcM{r^J) = e-4//3^ a, ■ a,. (2-1) 

Using the parameters listed in Table I, we can well reproduce masses of ordinary 
baryons and mesons as shown in Table II and III, and well reproduce or predict 
hadron properties, e.^., /Up ~ 2.75nm (exp. 2.78nm), — — l.SOnm (exp. — 1.91nm), 
/iA — — O.GOnm (exp. — O.Glnm), //^o ~ O.Slnm, /xq ~ — 1.84nm (exp. — 2.02nm) for 
the magnetic moment. Note here that the form of the linear-type confining potential 
V^onf is appropriate as an approximation of the Y-type linear three-quark potential 
indicated by lattice QCD,!^'!^ although the adopted value of the string tension a ~ 
0.56GeV/fm is smaller than the standard value a ~ 0.89GeV/fm. 



Table I. Parameters of the present constituent quark model with a linear-type confining potential. 



parameter ttIu, 


TTis string tension ct 


/3 


Wo 




adopted value 330MeV 


500MeV 0.56 GeV/fm 


0.5fm 


-572MeV 


240MeVxm^ 



Table II. Calculated masses of typical baryons obtained with the present constituent quark model. 



baryon 




N 


A 


A 


E 






N* 


calculated mass 


[MeV] 


939 


1235 


1064.8 


1129 


1324 


1539.5 


1462 


empirical mass 


[MeV] 


939 


1232 


1115.7 


1191 


1318 


1672.5 


1440 



Table III. Calculated masses of typical mesons obtained with the present constituent quark model. 



meson 




P 


K* 


K 


TT 


calculated mass [MeV] 


759 


759.3 


864.4 


457.9 


152 


empirical value [MeV] 


782.7 


775.5 


892 


496 


138 



Now, we perform an almost precise quark-model calculation with GEM for multi- 
quark systems and clarify whether each obtained state is a resonance or a continuum 
scattering state through the phase-shift analysis in the model calculation. For the 
five-quark system uudds with the energy E, the total wave function is given as 

^max 

^j.m{E) = ^fJSiE) + J2 b^j\E)'P^j}j,,{E,). (2-2) 
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Fig. 1. The five important sets of Jacobi coordinates among five quarks. Four ligfit quarks {u,d), 
labeled by particle 1 — 4, are to be antisymmetrized, while particle 5 stands for s quark. Sets 
c = 4, 5 explicitly contain two qq correlations, while sets c = 1 — 3 explicitly contain qq and qq 
correlations.'^ The NK scattering channel is treated with c = 1. 

The first term is the NK scattering component expressed by 

(E) = A,23,{ [ [<^f) (123)4^) (45)] . XL (Ri )] j.^,}. (2-3) 

The second term describes five-body degrees of freedom in the interaction-region 
amphtude which vanishes asymptotically. The amplitude is expanded using a nearly 
complete set of five-body eigenstates, = 1 — t'max)) constructed by 

diagonalizing the total Hamiltonian as {<Pj']j^j{Ey) \H\^^jJj^^{El)) = EySyyi. Here, 
each of <P^j}}^{Ey) is described as a superposition of L^-type five-body Gaussian 
basis functions written in all the Jacobi coordinates c = 1 — 5 as shown in Fig.l. By 
employing 15,000 five-body basis functions, i.e., z^max = 15,000, the eigenfunction 
set {^j^^jv/(^f )} forms a nearly complete set in the finite interaction region. 

The eigenfunctions (P^j}j^,j{E^) stand for discretized continuum states of the five- 
body system, and are called "pseudostates" in the scattering theory. It is known that 
most of the pseudostates do not actually represent resonance states but melt into 
non-resonant continuum states when the scattering boundary condition is imposed 
to the total wave function. 

We actually perform the precise quark-model calculation including all the pseu- 
dostate terms in (2.2), and calculate the phase shift 5 in the NK elastic scattering 
N+K^N-|-K for = |^,|^ channels, respectively. For J'^ = all the pseu- 
dostates located at — 450 MeV above the NK threshold melt into non-resonant 
continuum state when the NK scattering boundary condition is imposed, which 
means the strong coupling between those pseudostates and the NK scattering state. 
We show in Fig. 2 the calculated phase shift 6, and find no resonance for — 450 MeV 
above the NK threshold, i.e., 1.4 — 1.85 GeV in mass region around 0~''(154O). We 
find a sharp resonance for J'^ = ^ and a broad one for J'^ = ^ around 1.9 GeV in 
Fig. 2, although their energies are too high to be identified with the 0^(1540). For 
the J'^ = ^ states, we find no resonance up to 500 MeV above the NK threshold. 

Note that these results on the absence of low-lying pentaquark resonances with 
J'^ = ^ , I are consistent with the lattice QCD results in Refs.9) and 10). As for 
the presence of ^ penta-quark resonance around 1.9GeV, Ref.ll) shows a consistent 
lattice QCD result indicating a ^ penta-quark resonance around l.SGeV. 

To conclude, there is no five-quark resonance with J'^ = ^ , f below 1.85GeV. 
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Fig. 2. The calculated phase shift 5 for the five-quark system with J'^ — ^ (left) and — 
(right) in the NK scattering (N+K N+K) as the function of the energy measured from the 
NK threshold i5th ~ ttin + ttik- No resonance is seen around the energy of e+(1540). Instead, 
there is a sharp resonance state with J'^ = i and the mass of about l.QGeV. 



Instead, the quark-model calculation predict a five-quark resonance state of J'^ = | 
with the mass of about 1.9GeV and the width of i"' ~ 2.68MeV. 

§3. Four-body quark- model calculation for X°(3872) 

With GEM, we analyze four-quark systems ccqq for / = 0, 1 in a quark model 
to investigate X°(3872),'SJ which may have the tetraquark structure as ccuu\^ Note 
here that, if X'^(3872) is an / = 1 state, it is manifestly exotic and there should 
exist its isospin partners {ccud and ccdu) around 3.87GeV. For the calculation 
of four-quark charmed systems, we adopt the quark-quark interaction of Ref.13), 
which effectively includes the exchange effect of Nambu-Goldstone bosons. With 
rric — 1752MeV, rriu^d — 313MeV and ~ 555MeV, the quark model leads light 
hadron masses as nip ~ 772.8MeV, m^j ~ 696.3MeV and m,r — 148.7MeV, and the 
calculated (experimental) masses of charmed mesons are ttt-d — 1897.6(1867)MeV, 
mo* ^ 2017.1(2008)MeV, mj/^ ~ 3096.5(3097)MeV and m^^ ~ 2989.1(2980)MeV. 

For the four-quark calculation, we employ all the 18 sets of the Jacobi coordinates 
of the four-body system.^ We show in Fig. 4 the three important sets describing the 
cu and cu correlations (c = 1), the cc and uu ones (c = 2) and the cu and cu ones 
(c = 3). The four-body Gaussian basis functions are prepared in the Jacobi sets 

c = 1 — 18 and the four-body pseudostates 'l>^j} ^^^{Ey) are obtained by diagonalizing 
the four-quark Hamiltonian using nearly 13,000 Gaussian basis functions. For the 
phase shift calculation, we consider the scattering channels D'' -|- D*" and J/V' + w 
for I = 0, J'^ = 1+, and D° + D*° and J/V' + p iov I = I, J'^ = 1+ . 

For the / = states, we find a very sharp resonance, dominantly having the 
]2)Og*o component, slightly below the D^D**^ threshold. In this calculation, however, 
the u threshold energy is much lower than the experimental value by 86 MeV. 
We find that if the quark-quark interaction is artificially adjusted so as to reproduce 
the experimental value of the 3/ip to threshold, neither bound nor resonance state 
appears. In fact, the result seems rather sensitive to the quark-quark interaction, 
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and we need better interaction to obtain definite conclusions for / = states. 

For the 1=1 states, we obtain several pseudostates near the D'^D*'^ threshold 
region, but all of them disappear when the scattering boundary condition is switched 
on. The resultant scattering phase shift by the full coupled-channel calculation 
is given in Fig. 3, which shows no resonance behavior. We thus find no ccg^-type 
tetraquark resonance with 7 = 1 in mass region of 3.87— 4.0GeV. 

40 1 . ^ ^ ^ ^ r~, 




Fig. 3. (Left) Three important sets of the Jacobi coordinates in the scattering calculation on 
X°(3872). (Right) The phase shift in hadron-hadron scattering obtained by the four-body 
coupled-channel calculations for ccqq four-quark systems with / = 1, = 1^: the solid curve 
for J/ip + p-^ J/ip + p and the dotted curve for D° -|- D*° ^ D° -|- D*°. The scattering energy 
is measured from the J/i/j p threshold. The D^D*" threshold is at 45.3 MeV. No resonance is 
seen although a threshold cusp appears due to the couple-channel effect. 

In this way, the almost precise quark-model calculation using the Gaussian Ex- 
pansion Method (GEM) is a powerful tool to clarify the state properties of multi- 
quark systems, which is a theoretical search for exotic hadron resonances. 
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